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ABSTRACT 

We present our centimeter wavelength (1.4, 2.3 and 4.8 GHz) light curves of the afterglow of 
GRB 030329, which were obtained with the Westerbork Synthesis Radio Telescope. Modeling the 
data according to a coUimated afterglow results in a jet-break time of 10 days. This is in accordance 
with earlier results obtained at higher radio frequencies. However, with respect to the afterglow model, 
some additional flux at the lower frequencies is present when these light curves reach their maximum 
after 40-80 days. We show that this additional flux can be modeled with two or more components 
with progressively later jet breaks. From these results we infer that the jet is in fact a structured or a 
layered jet, where the ejecta with lower Lorentz factors produce additional flux that becomes visible 
at late times in the lowest frequency bands. We show that a transition to non-relativistic expansion of 
the fireball at late times can also account for the observed flux excess, except for the lowest frequency 
(1.4 GHz) data. 

Subject headings: gamma rays: bursts - radio continuum: general - radiation mechanisms: non- 
thermal 



1. INTRODUCTION 

Since the discovery of afterglow emission of gamma-ray 
bursts (GRBs) at X-ray, optical and radio wavelengths 
; llCosta et all Il997t iVan Paradii s et all Il997t iFrail et all 
, Il997j) . it has become clear that broad-band observations 
are needed to determine the physical processes produc- 
ing the afterglow emission in the context of the available 
, models, the most popular being the fireball m odel (e.g., 

■ iRees fc Mcszaros 1992; Mcszaro sfc R,eesill99lD . Obtain- 
ing the overall shape of the energy distribution and the 
time evolution of the GRB afterglow provides informa- 
tion about the intrinsic energy, both in electrons and in 

I magnetic fields, as well as about the matter into which 
, the G RB blasted its ejecta (see, e.g.. lWiiers fc Galamal 

■ Il999t) . Although optical observations alone can constrain 
I the value of some of these physical parameters, observa- 
tions covering the radio to X-ray wavelength regions are 
required to determine all of them. 

The self-absorption frequency, fa, of the after- 
glow broad-band spectrum can often be constrained 
by radio observations a t centimeter wavelengths (e.g., 
iWiiers fc Galamal I1999D . As the afterglow spectrum 
evolves, the two other characteristic break frequencies in 
its broad-band spectrum (the frequency at the peak flux, 
Vyn, and the cooling frequency, v^) enter the radio regime 
as well, although in practice the flux level at Vc is below 
the detection limit of current radio telescopes. The latter 
frequency can usually be determined from the available 
optical and X-ray data, which span the frequency range 
where the cooling frequency is found during the first days 
of the afterglow. These break frequencies and their time 
evolution uniquely determine the parameters that make 
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up the fireball model and its evolution. 

GRB 030329 is the closest gamma-ray burst discovered 
so far for which an af terglow has been fou nd^ . At a red- 
shift of z = 0.1685 l)Greiner et alJl2003t) . its afterglow 
was discovered in R = 12.4 m agnitude, just 67 minutes 
after the GRB itself l|Sato et al. 2003), about 100 times 
brighter than the average GRB afterglow. The bright- 
ness of the afterglow made it possible to study its evolu- 
tion for a long time and in detail over a broad range 
of frequencies, from X-ray to centimeter wavelengths. 
Furthermore, its proximity provided an excellent oppor- 
tunity to look for a supernova signature in both the 
light curve and the spectr um, as predicted by the coUap- 
sar model (Wooslcv 1993: MacFadven fc Wooslevl lT999). 
the currently favored progenitor model for long duration 
gamma-ray bursts. The resemblance between the super- 
nova spectrum distilled from the GRB 030329 afterglow 
and that of the energetic type Ic supernova SN1998bw 
(associated with GRB 980425, fCalama et al.l ll998') pro- 
vides strong support for the core colla pse of massive stars 
as th e cause for GRBs ijHiorth et al.ll2003l: iStanek et all 
l20?)l . 

Several authors have modeled the broad-band after- 
glow behavior with a standard fireball model for the 
afterglow. A first approximation shows excess flux 
(on top of the already bumpy light curve) after the 
first few days, most noticeable at the lower frequen- 
cies. JWillingalc et al. (2005) attribute the excess flux 
to the unde rl ying s upernova, bu t mos t authors (e.g. 
iBerger et all l|200l : iSheth et all l|2003|l : iTiengo et all 
|2(](]^) prefer a two-component jet model, where a slower 
jet is responsible for the extra emission appearing at op- 
tical wavelengths around 10 days after the burst. Even 
later time observations show a likely transition to the 
non-relativistic regime, estiinated a round 40 - 50 da ys 
after the burst ijTiengo et all l|200^ : iRajTeTaD l|20^ 

5 GRB 980425/SN 1998bw at z = 0.0085 was closer, but no af- 
terglow was found 
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Here we describe our radio monitoring campaign of 
this extraordinarily bright afterglow with the Wester- 
bork Synthesis Radio Telescope (WSRT) in the centime- 
ter waveband. In Section |21 we describe the data we ob- 
tained. In Section 13 we apply an afterglow model to the 
data, and in Section^lwe compare our results with radio 
data obtained by other groups. Finally, in Section |SJ we 
summarize our findings and draw our conclusions. 

2. DATA REDUCTION AND ANALYSIS 

Data were obtained with the WSRT, at 1.4, 2.3 and 
4.8 GHz. We used the Multi Frequency Front Ends ijTanI 
119911) in combination with the IVC-I-DZB backend** in 
continuum mode, with a bandwidth of 8x20 MHz. Gain 
and phase calibrations were performed with the calibra- 
tor 3C286, though sometimes 3C147 or 3C48 were used. 
Tablc^lists the log of the observations, all done in 2003. 
VLBI observations prevented us from using the WSRT in 
the second half of May, and observations were resumed 
in June, mostly at 4.8 GHz. At 2.3 and 1.4 GHz the 
observations suffered from confusion from nearby bright 
sources, causing the noise to be at least a factor of two 
above the theoretical limit. 

We checked our results for consistency by measuring 
the flux of several nearby point sources, which were as- 
sumed to be constant. In a few observations, we found 
these sources to be systematically dimmer, as indicated 
in the observation log; we therefore suspect that the flux 
derived in these observations for the afterglow is also be- 
low its real value. Although we could in principle scale 
these fluxes upward, we decided to ignore these observa- 
tions in our analysis, as the cause of these low flux levels 
is not clear. 

The 1.4, 2.3 and 4.8 GHz light curves are presented 
in Fig. n] The general trend of the light curves is that 
expected for the low-frequency part of a GRB afterglow: 
as long as the self-absorption frequency, u^, is higher than 
the observed frequency interval, the light curve rises since 
the frequency of the minimum electron injection energy 
Vrn moves toward lower frequencies. When both Vg, and 
i-',n pass the observed frequency interval (not necessarily 
at the same time), a turn-over in the light curve occurs 
and the flux falls off steeply. We have listed the precise 
temporal dependencies of fa, i^m: ^^c a-nd the peak flux 
Fai in Table 121 for a homogeneous circumburst medium 
as well as for a massive stellar wind, in which the external 
density p depends on distance r to the center as p oc r^^. 
The evolution in time of fa, fm and is plotted in Fig. 
12 Table 13 lists the dependencies of the flux on i/ and t 
(note the change in the spectral index when fm becomes 
less than fa)- 

3. APPLYING THE FIREBALL MODEL TO THE DATA 

We have modeled the data simultaneously in time and 
frequency, using a general broad-band afterglow model 
which includes a jet break and a transition to non- 
relativistic expansion of the fireball. The several power- 
law segments of the broad-band spectrum were connected 
smoothly in a way described in Appendix ^ Because 
-Fin, fa, fm and fc are functions of time, we need to quote 
them at some fixed moment, for which we choose the jet- 
break time. At earlier and later times the characteristic 



frequencies and the peak flux evolve in time according to 
Table El The resulting light curves are given in Table 13 
The transitions between the different regimes, marked by 
the jet-break time tj and the time tnr at which the fireball 
becomes non-relativistic, are treated as smoothly broken 
power laws as described in Appendix El 

Since our data show a large scatter, especially at early 
times due to scintillation, where they do not follow a 
smooth curve, we did not apply a, fit, but merely tried 
to obtain a best fit by eye. This ignores the scintillation, 
and it will also put some more emphasis on the 1.4 GHz 
light curve: a fit tends to follow the 4.8 GHz data 
points since they are more numerous, and will hence ig- 
nore the global trend seen in the 1.4 GHz light curve. 
The result of such an eye-ball broad-band fit is shown in 
Fig. n Note that the scintillation amplitude is largest 
in the first few days, after which it quickly declines (the 
large decrease in the 2.3 GHz light curve around 20 days 
could be an artifact in the data, since it does not show 
up in the other light curves). 

In this way, we obtained a value of 35 GHz for the 
electron injection frequency, fm, and a value of 13 GHz 
for the self absorption frequency, fa, at 10 days after the 
burst, which is the jet-break time. The flux at fm is 
about 61 mJy at that time. We find the electron index 
to be p = 2.1. The value of fc can only be determined 
with observations at higher frequencies; our dataset indi- 
cates that fc > 10^^ Hz , which is in agreement with the 
findings bv lSmith et alJ 1)20051) . who find that the rapid 
fall in their 350 GHz light curve can be attributed to the 
cooling frequency passing through their observing band. 
These results are so mewhat at odds with the findings by 
iBerger et al.l l)2003D . They obtain higher values for the 
characteristic frequencies fa (19 GHz) and fm (43 GHz), 
and the peak flux (96 mJy) at the jet-break time 
tj ~ 9.8 days. To investigate this further, we performed 
a fit which include s their data, as well as radio data from 
ISheth et al.lll2?ffll . This fit gives similar results to those 
obtained from our previous fit. 

From the obtained characteristic frequencies and the 
peak fiux we find an isotropic energy i^iso — 4.0 • 
lO^^f^^g erg, a circumburst density n ~ 0.56fc''i3 cm~'^, 
and the fractions of energy in the electrons Cc and mag- 
netic field eb of 0.25fyj3 and 0.49f "^3"' respectively. The 
cooling frequency can not be determined from the ra- 
dio observations at centimeter wavelengths, but we take 
" ^ 10^'^fc,i3 to compare our results with Bcrgcr et aO 



l)2003j) . The opening angle of the jet can be found to 
be 0j ^ 



0.38fc 13 rad (22°) from the jet-break time 



' See sect. 5.2 at|http : //www ■ astron ■ nl/wsrt/«srtGuide/node6 ■ html | 



of 10 days, which gives a beaming-corrected energy of 
Ecoi — 2.9 • lO^^f ~j^3^ erg. This energy is comparable to 
IBerger et al.l ()2003|) , the circumburst medium density we 
find is smaller, but Eo and eb are larger. 

At the turnover in the late-time light curves at 1.4, 
2.3 and 4.8 GHz we find an excess in flux compared to 
the model. We present two possible explanations for this 
behavior. The first one is a non-relativistic phase after 
im- ~ 80 days. TablcOshows that the light curves flatten 
in the transition to a non-relativistic phase. In Fig. ^o^is 
can see that this extension of the model fits the data at 
2.3 and 4.8 GHz well when one assumes a homogeneous 
xternal medium. However, the fiux at 1.4 GHz is over- 
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TABLE 1 

Log of the WSRT observations of GRB030329 in 2003 



observing dates 




integration time 


frequency flux 


error flux 




(days) 


(hours) 


(GHz) (mJy) 


(mJy) 



Mar 30.658 - 31.138 


1.414 


3.3 


1.4 


0.21 


0.07 


Mar 30.674 - 31.148 


1.427 


3.6 


2.3 


0.28 


0.05 


Mar 30.691 - 31.157 


1.440 


3.6 


4.8 


1.05 


0.03 


Mar 31.655 - Apr 1.110 


2.399 


3.6 


4.8 


5.98 


0.03 


Mar 31.672 - Apr 1.126 


2.415 


3.6 


2.3 


2.17 


0.05 


Mar 31.688 - Apr 1.142 


2.431 


3.3 


1.4 


0.63 


0.04 


Apr 2.650 - 3.088 


4.385 


4.0 


4.8 


3.64 


0.04 


Apr 2.694 - 3.132 


4.429 


4.0 


2.3 


0.79 


0.05 


Apr 2.738 - 3.149 


4.459 


3.3 


1.4 


0.43 


0.15 


Apr 4.644 - 5.083 


6.380 


4.0 


4.8 


4.89 


0.04 


Apr 4.688 - 5.127 


6.424 


4.0 


2.3 


1.21 


0.06 


Apr 4.732 - 5.143 


6.454 


3.3 


1.4 


0.37 


0.04 


Apr 5.641 - 5.818 


7.246 


4.2 


2.3 


1.00 


0.07 


Apr 7.072 - 7.138 


8.621 


1.6 


4.8 


4.21 


0.06 


Apr 7.636 - 8.135 


9.401 


12.0 


4.8 


3.96 


0.04 


Apr ll.DZo - ll.byo 


13.175 


1.6 


4.8 


6.42 


0.05 


Apr iz.bzz - iz.oUu 


14.230 


2.5 


4.8 


5.25 


0.05 


Apr iz.Doy - iz.o4o 


14.267 


2.5 


2.3 


0.26 


0.06 


Apr io.bUb - lo.9U5 


20.273 


7.0 


4.8 


6.96 


0.05 


Apr zU.bUU - zl.Uuo 


OO Q/1 /I 

ZZ.044 


Q Q 
O.O 


A Q 


O.OU 


U.U4 


A T^r 90 fi17 91 079 

/\pr zu.Dii - zi.u/z 


ZZ.ODU 


O.O 


9 

Z.O 


1 91 


U.UD 


Ar»r 90 fi'^'^ 91 087 


99 '^7^^ 


O.O 


1 4 




u.uo 


Ar^r 97 9S 0fi9 
/\pi z i .000 - ZO.UUZ 


9Q '^'^R 


o.y 




Q 1 7 




Ar^r 97 fiOfi 9S (X7^ 


9Q '?^7 


O.O 


9 

z.o 


2 32 




iViely Z.OUO - Z.O-l-U 






9 

z.o 


o.u^ 


u.uu 


ividy 0. i oz - o.ytjo 


o / .ouu 


O.O 




o.yu 






41.313 


0.7 


4.8 


8.05 


0.10 


Mav 9 81 3 - 9 841 


41.343 


0.7 


2.3 


3.11 


0.18 


Jun 9.714 - 9.963 


72.354 


2.5 


4.8 


4.59 


0.04 


Jun 9.737 - 9.944 


72.356 


2.5 


2.3 


3.84 


0.07 


Tun 1 6 567 - 1 6 645 


79.122 


1.9 


4.8 


4 3gb 


0.08 


Tun 1 7 7'^0 1 7 7Q9 


80.287 


1.0 


4.8 


3.79 


0.08 


Tun 1 S 777 1 SI Q 
J uii ±o. til ±o.oxy 


81.314 


1.0 


4.8 


3.74 


0.10 


Tiin 'in 4*^0 'iO 4.40 


92.951 


0.3 


4.8 


2.34 


0.18 


Till 1 fi'^9 1 74*^ 


94.214 


1.0 


4.8 


3.82 


0.10 


Jul z.4oo - Z.oz4 


oc^ n99 


n n 

u.y 


A Q 


9 ft 7b 


U.Uo 


Jul 3.548 - 3.587 


96.084 


1.6 


4.8 


3.57 


0.07 


Jul 5.892 - 5.933 


98.429 


1.0 


4.8 


3.09 


0.04 


Jul 19.580 - 19.854 


112.233 


6.6 


4.8 


2.27 


0.05 


Jul 23.673 - 23.843 


116.274 


4.1 


4.8 


2.43 


0.05 


Jul 29.519 - 29.644 


122.097 


3.0 


2.3 


3.00 


0.06 


Jul 29.675 - 29.800 


122.254 


3.0 


1.4 


1.72 


0.15 


Aug 2.546 - 2.595 


126.086 


1.2 


4.8 


2.13 


0.07 


Aug 2.619 - 2.668 


126.159 


1.2 


2.3 


0.80*= 


0.27 


Aug 2.692 - 2.740 


126.232 


1.2 


1.4 


1.93 


0.18 


Sep 13.535 - 13.701 


168.134 


4.0 


4.8 


1.28 


0.04 


Sep 16.196 - 16.250 


170.739 


1.3 


1.4 


2.29 


0.19 


Sep 28.161 - 28.285 


182.739 


3.0 


2.3 


1.95 


0.09 


Oct 11.291 - 11.458 


195.891 


4.0 


1.4 


1.40 


0.21 


Nov 29.248 - 29.491 


244.885 


5.8 


1.4 


1.15 


0.09 


Dec 1.128 - 1.288 


246.724 


3.8 


4.8 


0.85 


0.04 



''In days after the burst. The indicated time is the logarithmic average of the start 
and end of the integration. 

''The flux of surrounding point sources is consistently lower compared to other obser- 
vations. 



estimated at the latest times in this case (note that the 
symbols are larger than the error bars). The model with 
an external density gradient and a non-relativistic phase 
gives a better result at 1.4 GHz, but underestimates the 
flux at higher frequencies. 

The second explanation of the late-time flux excess is 
an extra component, which consists of an afterglow with 
a jet-break time later than 10 days. The resultant fit is 
also shown in Fig. ^ For the first component, which pro- 
duces the main flux at higher frequencies, the parameters 



are set as before except for ~ 48 mJy. The second 
component has tj ~ 30 days, i^m — 20 GHz, — 10 
GHz and Fm ~ 16 mJy at the jet-break time; so at 
t ~ 10 days the second component has ~ 35 GHz, 
I'll — 13 GHz and — 48 mJy. The electron index 
p 2.2 for both jets. The physical parameters we derive 
from these characteristic frequencies and peak flux are 

n ~ 0.82i/f_(3cm-3, Ce ~ 0.28i/^Yit and cb OASv'J^'^ 
for both components; for the first component we find 
Eiso ^ 4.0 • lO^V^^/gCrg, 0j ~ 0A2iy~ll^Tiid (24°) and 
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Fig. 1.— Eyeb all fits to the 4.8, 2.3 and 1.4 GHz WSRT data simultaneously, with 15 GHz VLA/ATCA/RT observations from 
IBerger et aTl 12003) for comparison. Note the large scatter around the model light curves during the first days. Left: The lines represent 
models with fm ~ 35 GHz, ~ 13 GHz and i^m — 61 mjy at tj = 10 days. The solid line corresponds to a model in which the fireball 
expands into a homogeneous medium and the non-relativistic phase of the fireball evolution starts after 80 days; the dotted line corresponds 
to the same model but without a non-relativistic phase, so it deviates from the solid line after 50 days; the dashed line corresponds to a 
model with a non-relativistic phase after 80 days and expansion of the fireball into a massive stellar wind. The peak frequency falls below 
the self absorption frequency at 17 days. From then on, the maximum of a light curve at a given wavelength marks the passing of the self 
absorption frequency. Right: A two-component fit to the data. The first component (dotted line, with a jet break time of 10 days) is 
responsible for the light curves until 50 days, while the second component (dashed line, with tj ~ 30 days) accounts for the later peak in 
the light curves. The combined light curve is shown as the solid line. 



TABLE 2 

The various temporal dependencies of the break frequencies and peak flux of the afterglow broad band spectrum. 
Before the jet-break time tj and after the non-relativistic timescale different scalings arise from a homogeneous 
circumburst medium or a stellar wind, when the external density p depends on distance r TO the center as p oc r~2. 
Between and tnr the external density profile does not influence the scalings. 







< Uc) i^a 


{Vm < Us, < Uc) 






Fm 


t < tj < tnr (homog 


eneous) 


t" 


t- 


(■ip+i)/i(p+i) 


(-3/2 


(-1/i 


t" 


t < tj < tnr (stellar 


wind) 


t-3/5 


t- 


3(p+2)/2(p+4) 


(-3/2 


(1/2 


(-1/2 






4 -1/5 


t- 


2(p+l)/(p+4) 


t-2 


tO 


t-1 


tj < tnr < t (homog 


eneous) 


(6/5 


t- 


(3p-2)/(p+4) 


t-3 


(-1/5 


(3/5 


tj < tnr < t (stellar 


wind) 


(-2/15 


t- 


(7p+6)/3(p+4) 


(-7/3 


tl 


(-1/3 
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TABLE 3 

The spectral and temporal flux dependencies in the different regimes of the broad-band afterglow spectrum. 





F {V< l^a) 


F (l/a 


<V <Vnx 


) F (^'nl 


< U < Uc) 




F (uc < u) 






t < tj < tnr (homogeneous) 


u'^-t^/-' 




^1/2 


i.-(p- 


l)/2 . j-3(p- 


i)/4 


u-pr^.t- 


(3p- 




t < tj < tnr (stellcir wind) 


u^-t^ 




tO 




l)/2 . j-(3p- 


l)/4 


y-P/'^-t- 


(3p- 


-2)/4 


tj <t < tnr 


U^-tO 


,.1/3 


t-1/3 


1/-(P- 


l)/2 .t-P 




v-P/^-t- 


P 




tj < tnr < t (homogeneous) 


u2.t-y'> 


iyl/3 


^8/5 




l)/2 . j-3(5p 


-7)/ 10 


v-P/^-t- 


(3p- 


-4)/2 


tj < tnr < t (stellar wind) 


!/2 .42/3 


!yl/3 


^4/9 


Jy-(P- 


l)/2 .t-(7p- 


5)/6 


v-P/^-t- 


(7p- 


-8)/6 




F <Vm 


) F {Un, 


< < l/a) F (^'m 


< U < Uc) 




F (uc < u) 






t < tj < tnr (homogeneous) 


r.2.tV2 




t5/4 




l)/2 . f-3(p- 


l)/4 


u-P/'^-t- 


C3p- 


-2)/4 


t < tj < tnr (stellar wind) 




1.5/2 


^7/4 




l)/2 . j-(3p- 


l)/4 


u-P/'^-t- 


(3p- 


-2)/4 


tj < t < tnr 


U^-tO 


^5/2 


tl 


J,-(p- 


l)/2 .f-p 




ij-P^- -t- 


P 




t) < tnr < t fhorriogcncous) 




,.5/2 


f 11/10 




l)/2 . ^-3(r,p 


-7)/10 


//-'"'- ■ / - 


Cip- 


-4)/2 


t] < ^nr < l- (sLc'llar wind) 






jil,(i 




l)/2.^ >7, 




w ■ ( 


[ 7/< 






1 10 100 1 10 100 

Time (days after burst) Time (days after burst) 



Fig. 2. The temporal evolution of the electron injection frequency, i/m, the self absorption frequency, i/a, and the peak flux, Fm- The 
upper panels show the evolution of Vm (solid line) and z^a (dash-dotted line), the lower panels the evolution of Fm. The left panels show 
a model in which the fireball expands into a homogeneous medium, with a jet-break time tj = 10 days, and the non-relativistic phase of 
the fireball evolution starts after tnr = 80 days. The right panels show a model in which the fireball expands into a massive stellar wind, 
with a jet-break time tj = 10 days, and the non-relativistic phase of the fireball evolution starts after tnr = 80 days. The dotted lines show 
where the breaks in the temporal behavior of the parameters occur at tj and tnr- 
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Ecor — 2.4 • 10^"i^(. 13 erg, while for the second compo- 
nent Eiso 9.1 • 10'5°i/^\(3 erg, 0j ~ O.rSiy^j^^ vad (42°) 

and E'cor ^ 2A-10^'^iy^ll^ erg. The data are well fitted by 
this two-component jet model. However, d ata at higher 
radio frequencies from iBerger et al.l l)2003D can not be 
fitted well in this model. 

4. DISCUSSION 

A similar procedure of fitting t wo component s with dif- 
ferent jet breaks was applied bv IBerger et alJ l|200 3) to 
explain the break in the early-time optical (and X-ray) 
light curve. The underlying mechanism involves two jet- 
ted outflows, one with a small opening angle and a high 
Lorentz factor that produces the early-time light curve 
(with tj ~ 0.5 days), and a jet with a larger opening 
angle and lower Lorentz factor that carries the bulk of 
the energy and produces the later-time light curve (with 
tj ~ 10 days). The WSRT observations at 2.4 days af- 
ter the burst have values for the flux that are well above 
the theoretical curves (see Fig. QJ. We investigated the 
possibility that these are signatures of the jet that pro- 
duces the early-time optical light curve. However, with 
the constraints on the parameters from the optical and 
X-ray observations, it is not possible to fit the early ra- 
dio observations with this jet with a jet-break time of 

0. 5 days. A better explanation for these observations is 
scintillation. 

From the result of our two-component model fit, we 
can conclude that, besides the jets with tj ~ 0.5 and 
tj ~ 10, another jet is present with even larger opening 
angle, that powers the late-time {t > 50 days) light curve 
and is therefore best visible at the very low frequencies 
observed here. However, it may be that the total jet 
(which possibly includ es the fir st narrow jet as well) is 
structured (e.g. Meszaros et al.l ll998: Rossi et al. 20021 
and that the Lorentz factor T decreases toward the edge 
of the jet-cone. Alternatively, the outflow consists of a 
layered jet, where shocks with lower F follow the faster 
ones as they run into the surrounding medium. In both 
cases, one expects that the low Lorentz factors dominate 
at low frequencies and late times, and that the jet break 
occurs later at progressively lower frequencies. 

The multiple component model is certainly not satis- 
factory: it does not give a good fit to the data at radio 
frequencies above 4.8 GHz at late times. Our model in 
which a transition to a non-relativistic phase of the fire- 
ball occurs after 80 days gives a better broadband radio 
fit except for the data at 1.4 GHz. This transition to 
a non-relativistic phase is also seen in VLA radio obser- 
vatio ns by Frail et al. ( 2005) and at X-ray frequencies 
bv iTiengQ et alJ l(2004ll . although their estimate for the 
time at which this transition occurs is lower than ours, 

1. e. ~ 50 and 44 days respectively. Our low frequency 
radio data can not be fitted well by applying this low 
value for the non-relativistic transition. 

Although the model in which a transition to a non- 



relativistic phase of the fireball occurs, gives the best 
broadband radio fit, the value oitar we find is much lower 
than theoretical estimates done bv lGranot et al.l 1)20051) 
and lOren et al.l l)2004|) . based on determin ations of the 
evolution of the image size of GRB 030329 fT avlor et a,lJ 
(|2005l) ). We estimate t„r ~ 209{r^^y/^ ~ 
168^^. ]^3^ days, which is a factor of 2 higher than the 80 
days we get from modeling the centimeter light curves. 
This discrepancy could be solved by fitting the broad- 
band afterglow light curves of GRB 030329 simultane- 
ously with the evolution of its image size. 

5. SUMMARY AND CONCLUSIONS 

Our data confirm the picture of a second jet in the 
afterglow of GRB 030329, that manifests itself around 
tj — 10 days. However, the flux level around the time 
when the low frequency light curves peak is higher than 
tha t predicted by the two-component afterglow model 
(cf. IBerger et al]l200l . Adding a third component with 
a later tj , we can account for this excess flux. Taking into 
account the early jet break, seen most clearly at optical 
wavelengths, we suggest that one is actually seeing the 
resul t of several blast waves with a range in Lore ntz fac- 
tors ijRees fc Meszaros..l998t iGranot et al.ll2003)) , some- 
thing whic h comes quite naturally in the coUapsar model 
for GRBs llMacFadven et alJ I200H iRamirez-Ruiz et alJ 
2002HZhang et aljl2003D . aiid was already suggested by 
Sheth et al.l 1) 2003). However, their high frequency data 
was unable to distinguish jet breaks at later times. Our 
later time low frequency radio data show such a late- 
time jet break, corresponding to a lower Lorentz factor, 
and therefore support a layered or structured jet for the 
afterglow of GRB 030329. 

An alternative explanation is the transition to a non- 
relativistic phase of the fireball. This model gives a good 
fit to the data at 2.3 and 4.8 GHz, but overestimates 
the flux at 1.4 GHz. This overestimation can be caused 
by the method of smoothly broken power laws as de- 
scribed in Appendix ^ We did not take into account 
the jet that is pointing away from us, and this can pos- 
sibly give extra flux at late times when the jet becomes 
non-relativistic and spherical. Continuation of observa- 
tions at late times at low radio frequencies and more 
detailed physical models can diagnose the cause of this 
discrepancy more closely. 
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thors acknowledge benefits from collaboration within the 
Research Training Network "Gamma- Ray Bursts: An 
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HPRN-CT-2002-00294. 



APPENDIX 

BROADBAND SPECTRUM AND LIGHT CURVE MODELING 

The dominating radiation mechanism for GRB afterglows is synchrotron emission. The broadband synchrotron 
spectrum is determined by the peak fiux and three break frequencies, namely the synchrotron self- absorption frequency, 
i^a, the frequency that corresponds to the minimal energy in the electron energy distribution, i^m, and the cooling 
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frequency, Vc, that corresponds to electrons that lose their energy quickly by radiation. The time evolution of these 
four parameters gives the evolution of the spectrum and thus light curves at all observing frequencies. 

The relativistic electrons emitting the synchrotron radiation are accelerated at the shock front to a power-law 
distribution of energies, 7V(7e)(i7e oc 7~^(i7e for 7e > 7m. For electrons that cool on a timescale smaller than the 
dynamical timescale, the energy distribution is steeper, N{^e)dle oc 7~^~^d7e- The broadband energy spectrum is 
found by connecting these two energy distributions at the cooling Lorentz factor 7c, and then integrating the single- 
electron synchrotron spectrum over the distribution function. Synchrotron self-absorption is taken into account by 
calculating the flux as follows: 

ii_(^)-i[l_exp(-^)], (Al) 



with the emission coeff icient, au the absorption coefficient, and D the distance. A detailed description of our 
modeling is presented in iVan der Horst et alJ l)2005j) . 

The evolution of the characteristic fr equencies and the peak flux is given in Table El We assume that after the 
jet-break time the jet spreads sideways (|Rhoadsl l)199fll) ). until it becomes spherical approximately at the same time 
the fireball becomes non-relativistic and enters the Sedov-Von Neumann- Taylor phase of the evolution. The transitions 
between the different regimes, marked by the jet-break time tj and the time inr at which the fireball becomes non- 
relativistic, are treated as smoothly broken power laws. We introduce a smoothening parameter s and take v.^ and 
as examples: 



3/2-s 



2-s 



-1/2-s 



2-s 



1/5-sN 



3 s 



1/s 



-1/s 



(A2) 



(A3) 



Expressions for and F^-a are similar. 

We choose to have a smoothening parameter of s = 5 for every transition. However, these transitions will probably be 
different and this can only be accounted for in detailed hydrodynamical modeling of the fireball. So this smoothening 
parameter gives a an uncertainty which may account for the discrepancy seen at 1.4 GHz in the case of one jet with a 
transition to the non-relativistic phase at inr — 80 days. 
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